Temporal evolution of Seebeck coefficient in an ac driven strongly correlated quantum 
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We study the response of the thermopower of a quantum dot in the Kondo regime to sinusoidal 
displacement of the dot energy level via a gate voltage using time dependent non-crossing approxi- 
mation and linear response Onsager relations. Instantaneous thermopower begins to exhibit complex 
fluctuations when the driving amplitude is increased at constant driving frequency. We also find that 
the time averaged thermopower decreases steadily until it saturates at constant driving amplitude 
as a function of inverse driving frequency. On the other hand, time averaged thermopower is found 
to be quite sensitive to ambient temperature at all driving frequencies for large driving amplitudes. 
We discuss the underlying microscopic mechanism for these peculiarities based on the behaviour of 
the dot density of states. 

PACS numbers: 72.15.Qm, 85.35.-p, 71. 15. Mb 



I. INTRODUCTION 

Single electron transistors consisting of a nanodevice 
placed between electrodes have received much attention 
recently since they are believed to carry great promise as 
a future replacement for conventional MOSFET which 
will reach its limit at 10 nm in 2016^. Unparalleled pre- 
cision enabled by nanotechnology revolution brought im- 
mense momentum to this field. Determination of the 
switching behaviour of these devices in the presence of 
strong electron-electron correlations is a requirement for 
their utilization in electrical circuits. Moreover, investi- 
gation of transient electron dynamics^ should also help 
fine tuning of quantum computerSiS, as well as single elec- 
tron gunsi. 

Initial investigations on the effect of Kondo correla- 
tions on the transient current flowing immediately after 
after abrupt changing of the gate or bias voltage^"— un- 
covered different time scales*'^. Since it takes Kondo res- 
onance considerable amount of time to get fully formed-, 
strong interference between the Kondo resonance and the 
Van Hove singularities in the contacts' density of states 
results in damped oscillations in the long timescale for 
an asymmetrically coupled systemic. Better agreement 
between the experiments and theoretical predictions can 
be attained by using realistic density of states for the 
contacts with the help of ab initio methods^^'^^. 

Transport experiments can derive additional benefit 
from the measurement of thermopower (Seebeck coeffi- 
cient) S since its sign is a benchmark providing direct 
information about the alignment of the impurity orbitals 
with respect to the Fermi level of the contacts. Conse- 
quently, this provides unambiguous evidence about the 
existence of Kondo resonance at the Fermi level of the 
contacts. Recent experiments measuring thermoelectric- 



ity in molecular junctions yielded key signatures confirm- 
ing these predictions^^"— . 

Previous theoretical studies on this set-up mainly fo- 
cused on steady-state behaviour of thermopower. Dif- 
ferent groups independently observed that the position 
of the dot level is critical to determine the sign of ther- 
mopower in Kondo regimei^ii^. In the first attempt to 
study the transient behaviour of thermopower™, it was 
uncovered that the inverse of the saturated decay time 
of thermopower to its steady state value is equal to the 
Kondo temperature when the dot level is abruptly moved 
to a position close to the Fermi level such that the Kondo 
resonance is present in the final position. 



It is quite interesting to investigate the effect of chang- 
ing the left and right tunnel couplings of the dot to 
the electrodes sinusoidally with a phase lagSS. If the 
phase lag is nonzero, this results in asymmetrical cou- 
pling. Another interesting situation arises when the dot 
energy level is subjected to time dependent ac pertur- 
bations in symmetrical coupling. This enables to moni- 
tor the evolution of thermopower in real time when the 
Kondo temperature changes constantly as a result of the 
oscillatory motion of the dot energy level between the 
leads^^-^^. This type Kondo shuttling has been made 
possible with experimental advances on controlling dot- 
lead coupling2^i2i. We previously studied the behaviour 
of current in this set-up and found that the time av- 
eraged conductance exhibits considerable deviation from 
its monotonic decrease as a function of applied bias when 
the bias is equal to the driving frequency^. It is our aim 
in this paper to study both the instantaneous and time 
averaged values of the thermopower when the dot level 
is driven sinusoidally by means of a gate voltage. 
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II. THEORY 

The physical set-up of a single electron transistor con- 
sists of a single discrete degenerate energy level Sdot cou- 
pled to continuum electrons in the electrodes. This is a 
quantum impurity problem and it can be represented by 
the Anderson Hamiltonian 

H{t) = ^ £dot{t)na + ^ £knkaa + \^^U^^a>n„n„, + 

(7 kaa 

^ [K<(efc«,OcL.c,-f H.c] , (1) 



adequately, where cj. (co-) and 



'kaa 



(ckaa) with a=L,R 



create (annihilate) an electron of spin a in the dot en- 
ergy level and in the left(L) and right (R) electrodes re- 
spectively. The ria- and nkaa are the number operators 
for the dot level and the electrode a. Va are the tunnel- 
ing amplitudes between the electrode a and the quantum 
dot. The Coulomb repulsion energy U is assumed to be 
infinity hence the occupation of the dot level is strictly 
restricted to unity. In this paper, we will adopt atomic 
units where h = kQ = e = 1. 

In our procedure, we perform slave boson transforma- 
tion for the Anderson Hamiltonian. The electron opera- 
tors on the dot are rewritten in terms of a massless (slave) 
boson operator and a pseudofermion operator as 



subject to the requirement that 



(2) 



(3) 



The last requirement is essentially a bookkeeping tool 
that prevents the double occupancy in the original Hamil- 
tonian when [/ — > cxD. The main advantage of using this 
transformation is that we can now safely drop the quar- 
tic Hubbard term since its effects are now incorporated in 
other terms anyway. The transformed Hamiltonian then 
becomes 



[^krikaa +Va{eka,t)cl^^b^ fa +}i.c. ,(4) 

kaa- 

where f^{fa) and b'^°-3{h) with Q!=L,R create(annihilate) 
an electron of spin a and a slave boson on the dot re- 
spectively. 

Assuming no explicit time dependency for the hopping 
matrix elements, the coupling of the quantum dot to the 
electrodes for a symmetrically coupled system can be ex- 
pressed as r(e) — f/9(e), where f is a constant given by 
f = 27r|y(e/)p and p{e) is the density of states function 
of the electrodes. 

The retarded Green function can now be rewritten 
in terms of the slave boson and pseudofermion Green 



functional^ as 
G^(t,ii) 



+ G^seudoit^h)B''{h,t)] (5) 



-9r 




'dot 




F„-F„,cos(nO 




-9r 



FIG. 1: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 

It is a quite a cumbersome and highly nontrivial task 
to compute the values of these double time Green func- 
tions in real time as the dot level is in constant mo- 
tion. We accomplish this by solving coupled integro- 
differential Dyson equations in a two-dimensional carte- 
sian grid. However, the self-energy of the pseudofermion 
and slave boson must be built into these equations in 
order to be able to obtain a solution. To this end, 
we invoke non-crossing approximation(NCA) to express 
the pseudofermion and slave boson self-energies^i^. As 
the name implies, NCA ignores higher order corrections 
where propagators cross each other. Consequently, it 
provides erroneous results in finite magnetic fields and 
when the ambient temperature is an order of magnitude 
smaller than Kondo energy scale. We will avoid these 
cases in this paper. Main advantage of NCA is that it is 
known to give accurate results for dynamical quantities. 
Upon solving the Dyson equations properly, their values 
are stored in a square matrix which is translated one step 
at a time in diagonal direction. 

In linear response, the conductance of the device is 
given by 



G{t) 



Lujt) 
T ' 



(6) 



and the Seebeck coefficient or the thermopower is ex- 
pressed as 



S{t) 



(7) 
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where the Onsager coefficients are 

/,„(/_^,./|i»-<-..)r(.)c.(M.)^tl)(8) 




-0.15 



and 

z„(/*j,,/ie'.<.-'..rWG-(M.)^fl)<s, 

where /(e) = ^^^/^t is the Fermi-Dirac distribution with 

Instantaneous conductance exhibits oscillations be- 
tween zero occurring at the lowest position of the dot 
level and a maximum taking place at its closest location 
to the Fermi level. We previously analyzed the conduc- 
tance for this set-up in detail^^ and presented the time- 
averaged density of states for the dot as well. 
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FIG. 2: This figure shows the instantaneous thermopower 
S{t) immediately after the gate voltage has been turned on 
for driving amplitudes of 2. OF (blue dashed), 2.5F (red solid) 
and 3.0r (black dot dashed) at r=0.003F and r2=0.14F in 
linear response. 



III. RESULTS AND DISCUSSION 

In this work, we will investigate the behaviour of 
Eq. ([7]) for a quantum dot in Kondo regime whose en- 
ergy level is driven sinusoidally via a gate voltage. We 
will report our results in atomic units therefore one needs 
to multiply our thermopower values with ks/e to obtain 
SI value in terms of Volts/Kelvin. This conversion fac- 
tor turns out to be 0,86x10"'*. This sinusoidal motion 
results in periodic modulations of the main Kondo res- 
onance and its satellites. Time averaged thermopower 
has been previously studied in this system neglecting 
the Kondo effectS,. We would like to make this model 
more realistic by including the Kondo effect in this anal- 
ysis since strong electron correlation effects play a crucial 
role in confined nanostructures like quantum dots. We 
will only consider the linear response instantaneous ther- 
mopower due to the validity of Onsager relations in this 
regime. 

Kondo effect has long been a trademark of many body 
physics taking place at sufficiently low temperatures due 
to a hybridization between the net spin localized inside an 
impurity and the continuum electrons of a nearby metal. 
Its main manifestation is a sharp resonance around the 
Fermi level of the metal arising when the impurity level 
is located below the Fermi level. The linewidth of this 
Kondo resonance is roughly equal to a low energy scale 
called the Kondo temperature which is denoted by Tk 
and expressed as 

In Eq. (|10p D is the half bandwidth of the conduction 
electrons while F = rp(e/). 



We will first consider the behaviour of the instanta- 
neous value of the thermopower when the dot energy level 
is modulated sinusoidally via a gate voltage. We can rep- 
resent the instantaneous behaviour of the dot level as 

edot{t)^-5r-Acos{nt). (11) 

where A is the driving amplitude and 51 is the driving 
frequency. The dc value of the dot energy level was cho- 
sen as — 5F to ensure that the dot level lies sufficiently 
below the Fermi level such that the Kondo temperature 
is much lower than the ambient temperature. In princi- 
ple, it is possible to choose a lower dot energy level as 
dc value. However, it would be problematic with our 
parameters because we would need to increase the driv- 
ing amplitudes even higher in that case to approach the 
Fermi level sufficiently and the dot energy level would 
go beyond the bandwidth which is at — 9F in its lowest 
position. This type of modulation induces a continuous 
change in Kondo temperature Tk since its value depends 
on the dot energy level in Eq. pU)) . Consequently, the 
shape of the Kondo resonance pinned to the Fermi level of 
the contacts constantly changes. The parabolic structure 
of the density of states of the contacts and the temper- 
ature gradient applied to them is depicted schematically 
in Fig. [T] alongside with the application of gate voltage 
to the dot energy level. 

Instantaneous thermopower is shown in Fig.[2]for three 
different driving amplitudes A immediately after the gate 
voltage has been turned on at constant driving frequency 
and ambient temperature. Since the dot level is oscillat- 
ing sinusoidally, the value of the Kondo temperature is 
fluctuating alongside with it between zero and a maxi- 
mum value. These maximum values are Tk = 0, 0022r, 
Tk = 0,00055F and Tk = 0,0001 IF for driving ampli- 
tudes of 3F, 2, 5F and 2F respectively. In this figure, 
thermopower slowly acquires a positive value as the dot 



4 



level approaches the Fermi level for all driving ampli- 
tudes. For the smallest driving amplitude, it reaches a 
maximum positive value at its closest point to the Fermi 
level before it starts decreasing again as the dot level 
begins to go down. 

However, this smooth behaviour changes abruptly 
when we start increasing the driving amplitude. When 
the driving amplitude becomes 2.5r, value of the instan- 
taneous thermopower initially acquires a positive value 
as in the previous case, however it suddenly starts to de- 
crease as the dot level approaches the Fermi level and 
dips into negative territory before it rebounds and in- 
creases again as the dot level gets away from the Fermi 
level. When the driving amplitude is increased even fur- 
ther, instantaneous thermopower starts decreasing earlier 
and more importantly, it dips deeper into negative terri- 
tory before it rebounds as the dot level is pushed back. 
This fluctuating behaviour is repeated over many cycles 
since the dot level's motion follows the gate voltage and 
is strictly sinusoidal consequently. 
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FIG. 3: This figure shows tlie instantaneous thermopower 
S{t) immediately after the gate voltage has been turned on 
for driving amplitudes of 2. OF (blue dashed), 2.5r (red solid) 
and 3.0r (black dot dashed) at T^O.OOSr and Sl^O.lOF in 
linear response. 

Fig. [3] shows the instantaneous thermopower for the 
same parameters used in Fig. [21 but the gate voltage is 
applied with a smaller frequency this time. An obvious 
consequence of this on the instantaneous thermopower is 
the longer oscillation periods for all driving amplitudes. 
Besides this trivial effect, we notice that the smooth os- 
cillation behaviour for the smallest driving amplitude 
did not change much. On the other hand, there are 
significant changes for the other two larger driving am- 
plitude cases compared to the larger driving frequency. 
Their fluctuation behaviour seems qualitatively similar 
but they both dive much deeper into negative territory 
as the dot level comes close to the Fermi level. This effect 
is quite pronounced for both cases. 

Having observed these peculiarities, one would be nat- 
urally inclined to find out whether these conclusions are 
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FIG. 4: This figure shows the value of the thermopower av- 
eraged over a period as a function of inverse driving frequency 
at r=0.003r (black dashed) and r=0.002r (blue dot dashed) 
for driving amplitude of 2.5r. 



valid for other driving frequencies and ambient temper- 
atures. In order to pinpoint this issue, we calculated 
the time average value of the instantaneous thermopower 
over a full period for two different ambient temperatures 
at several different driving frequencies. Time averaged 
thermopower for the smallest driving amplitude was not 
sensitive to the changes in temperature so we are not re- 
porting it in the following. We will elucidate its reason 
later. 

Fig. m depicts the time averaged value of the instan- 
taneous thermopower as a function of inverse driving 
frequency for two different ambient temperatures at a 
constant driving amplitude of 2.5F. The first conclusion 
that can be drawn from this figure is that the absolute 
value of the time averaged thermopower is always dif- 
ferent at all driving frequencies when the ambient tem- 
perature changes. Moreover, absolute value of the time 
averaged thermopower starts to increase for both ambi- 
ent temperatures as the driving frequency is decreased 
but it saturates at a certain driving frequency and stays 
constant thereafter. 

This is a quite peculiar and interesting observation and 
in order to test whether it is valid for other driving am- 
plitudes as well, we performed the same analysis for the 
largest driving amplitude of 3. OF keeping all the other pa- 
rameters constant. We do not increase the driving ampli- 
tude beyond this value to avoid entering mixed valence 
regime where the Kondo and Breit-Wigner resonances 
overlap. We are only interested in investigating the effect 
of Kondo resonance on thermopower in this paper. The 
result for this driving amplitude can be seen in Fig. [S] 
It is clear from this figure that the conclusion we had 
drawn is indeed valid here too. Namely, time averaged 
thermopower is quite sensitive to ambient temperature 
and it saturates when the driving frequency is decreased. 
However, there is an important difference. Absolute val- 
ues of the time averaged thermopower at both ambient 
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temperatures and all driving frequencies are greater than 
the previous case in Fig. 3] where the driving amplitude 
was 2. sr. This is also an important result which we will 
dwell on below. 

We now want to provide an unambiguous microscopic 
explanation to the results we presented above. It is clear 
that the development of the Kondo resonance plays a 
crucial role in transport of this prototypical device. Evo- 
lution of the Kondo resonance is incorporated in the dot 
density of states and the relation between the dot den- 
sity of states and the thermopower at low temperatures is 
given by the Sommerfeld expansion. It can be expressed 
as 



S{T) = - 



TT^T dA. 
3A{0,T)~de 



|£=0 



(12) 



in atomic units. In this equation, A{0,T) is the value of 
the dot density of states at Fermi level of the contacts 
and ^ is its derivative. It has been shown that the for- 
mation of the Kondo resonance occurs slightly above the 
Fermi level^^ making the derivative of the dot density of 
states at Fermi level positive at ambient temperatures 
close to or less than the Kondo temperature. This will 
give rise to a negative thermopower in this regime. On 
the other hand, Breit-Wigner resonance around the dot 
energy level is the only prominent feature of the dot den- 
sity of states in the absence of Kondo resonance at ambi- 
ent temperatures much greater than the Kondo temper- 
ature. Consequently, the thermopower would be positive 
if the dot level is positioned sufficiently away from the 
Fermi level since the dot density of states at Fermi level 
will lie in the long tail of the Breit-Wigner resonance. 
Thus, a shift in the sign of the thermopower by tuning 
the temperature is a signature of the emergence of the 
Kondo resonance. 
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FIG. 5: This figure shows the value of the thermopower av- 
eraged over a period as a function of inverse driving frequency 
at r=0.003r (black dashed) and r=0.002r (blue dot dashed) 
for driving amplitude of 3. OF. 

When the driving amplitude is 2. OF, the dot level can- 
not approach the Fermi level to be able to develop the 



Kondo resonance. Therefore, T^ Tk for all dot levels in- 
volved in this oscillation and the thermopower never dips 
into negative territory. It starts from zero when the dot 
level is farthest from the Fermi level because the den- 
sity of states is essentially zero everywhere around the 
Fermi level. However, it gradually acquires a positive 
value as the dot level approaches the Fermi level because 
the tail of the Breit-Wigner resonance starts to reach the 
Fermi level causing a negative slope there. Thermopower 
reaches a maximum positive value before the dot level 
starts to retreat. The values are not sensitive to temper- 
ature because Kondo resonance cannot be formed at all 
at any stage of this motion. 

The oscillation of the dot energy level covers a larger 
energy range as the driving amplitude is ramped up to 
2. sr. The value of the thermopower initially follows the 
value obtained when the driving amplitude was 2. OF but 
it starts to deviate significantly when the dot level ap- 
proaches the Fermi level because the Kondo resonance 
begins to emerge there. That is accompanied with a 
change in slope of the density of states of the dot and 
consequently a change in sign of the thermopower. The 
thermopower sees its minimum value when the dot level is 
at its closest point to the Fermi level because the Kondo 
resonance is at its most developed form here. As the 
dot level starts to retreat, the Kondo resonance is gradu- 
ally suppressed hence the thermopower becomes positive 
after a while. 

Similar conclusions hold for the largest driving ampli- 
tude of 3. OF except for the fact that the instantaneous 
thermopower dips lower into negative territory since the 
dot energy level manages to get even closer to the Fermi 
level hence the Kondo resonance is more developed. 

The saturation behaviour of the time averaged ther- 
mopower at low driving frequencies is also closely re- 
lated to the development of the Kondo resonance. As 
the driving frequency is lowered, dot energy level begins 
to spend more time near the Fermi level. This in turn 
gives the Kondo resonance more time to develop since its 
full formation takes considerable amount of time^. Con- 
sequently, time averaged thermopower values start to in- 
crease in negative territory. However, when the driving 
frequency is lowered below a certain value, we observe 
that the time averaged thermopower saturates. This is 
simply because the dot level is allocated sufficient time 
in the vicinity of Fermi level anymore. 

Dependence of the time averaged thermopower curves 
to ambient temperature is quite interesting too and re- 
quires further elaboration. It is clear from both Fig. |4] 
and Fig. [S] that the difference between the time averaged 
values obtained at the same driving frequency increases 
as the driving frequency is reduced until it stays constant 
below a certain frequency. This behaviour is further tes- 
timony to the crucial role development of the Kondo reso- 
nance plays in thermopower. At elevated driving frequen- 
cies, time averaged thermopower is not sensitive to ambi- 
ent temperature because the dot level cannot spend suffi- 
cient time near the Fermi level. Naturally, the Kondo res- 
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onance cannot be formed adequately hence a lack of sensi- 
tivity to ambient temperature. As the driving frequency 
is reduced, the difference becomes more pronounced and 
stays constant after the Kondo resonance gets the chance 
to develop itself fully below a certain frequency. 

IV. CONCLUSION 

In conclusion, we studied in linear response both the 
instantaneous and time averaged values of thermopower 
for a quantum dot sinusoidally driven by a gate voltage. 
We found that the instantaneous thermopower exhibits 
complex fluctuations as the driving amplitude is ramped 
up. We interpreted the origin of dipping into negative 
territory as a sign of development of the Kondo reso- 
nance at Fermi level. We studied the behaviour of the 
time averaged thermopower as a function of inverse driv- 



ing frequency and observed that it tends to saturate at 
low driving frequencies. Moreover, the results depend on 
ambient temperature quite sensitively. We believe that 
the scenario we discussed in this paper has great poten- 
tial to be realized experimentally. The temperature gra- 
dient between the contacts can be achieved with laser 
irradiation and the sinusoidal motion of the dot level can 
easily be realized with an ac gate voltage. We also think 
that state of the art ultrafast pump-probe technique o^^i"^" 
could make recording thermopower in real time possible. 
Therefore, we hope to create motivation and experimen- 
tal progress in this field with our novel results. 
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